Sex and statin-related genetic associations at the PCSK9 gene locus: results of genome-wide association meta-analysis.

BACKGROUND: Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a key player of lipid metabolism with higher plasma levels in women throughout their life. Statin treatment affects PCSK9 levels also showing evidence of sex-differential effects. It remains unclear whether these differences can be explained by genetics. METHODS: We performed genome-wide association meta-analyses (GWAS) of PCSK9 levels stratified for sex and statin treatment in six independent studies of Europeans (8936 women/11,080 men respectively 14,825 statin-free/5191 statin-treated individuals). Loci associated in one of the strata were tested for statin- and sex-interactions considering all independent signals per locus. Independent variants at the PCSK9 gene locus were then used in a stratified Mendelian Randomization analysis (cis-MR) of PCSK9 effects on low-density lipoprotein cholesterol (LDL-C) levels to detect differences of causal effects between the subgroups. RESULTS: We identified 11 loci associated with PCSK9 in at least one stratified subgroup (p < 1.0 × 10-6), including the PCSK9 gene locus and five other lipid loci: APOB, TM6SF2, FADS1/FADS2, JMJD1C, and HP/HPR. The interaction analysis revealed eight loci with sex- and/or statin-interactions. At the PCSK9 gene locus, there were four independent signals, one with a significant sex-interaction showing stronger effects in men (rs693668). Regarding statin treatment, there were two significant interactions in PCSK9 missense mutations: rs11591147 had stronger effects in statin-free individuals, and rs11583680 had stronger effects in statin-treated individuals. Besides replicating known loci, we detected two novel genome-wide significant associations: one for statin-treated individuals at 6q11.1 (within KHDRBS2) and one for males at 12q24.22 (near KSR2/NOS1), both with significant interactions. In the MR of PCSK9 on LDL-C, we observed significant causal estimates within all subgroups, but significantly stronger causal effects in statin-free subjects compared to statin-treated individuals. CONCLUSIONS: We performed the first double-stratified GWAS of PCSK9 levels and identified multiple biologically plausible loci with genetic interaction effects. Our results indicate that the observed sexual dimorphism of PCSK9 and its statin-related interactions have a genetic basis. Significant differences in the causal relationship between PCSK9 and LDL-C suggest sex-specific dosages of PCSK9 inhibitors.

The protein "proprotein convertase subtilisin/kexin type 9" (PCSK9) regulates the levels of low-density lipoprotein cholesterol (LDL-C) in blood, and thus, contributes to the risk of cardio-vascular diseases. Women tend to have higher PCSK9 plasma levels throughout their life, although the difference is smaller in patients under LDL-C lowering medication (e.g., statins). We investigated the interplay of genetics, statin-treatment and sex, using combined data from six European studies. We detected 11 genetic regions associated with PCSK9 levels, of which one was specific for women (at SLCO1B3, a statin-transporter gene), and three were specific for men (e.g., ALOX5, encoding a protein linked to chronic inflammatory diseases such as atherosclerosis). We also tested if statin use changed the genetic effect and found five genes only associated with PCSK9 levels in untreated participants. Variants in the gene encoding PCSK9 were most strongly associated and had heterogeneous effects in dependence on statin treatment and sex: On one hand, there were genetic variants with stronger effects in men than women. Those variants are also linked to sex-differential gene expression of PCSK9. On the other hand, there were also variants with treatment-depending effects, linked to protein structure and functionality of PCSK9. This indicates that the observed sexual and treatment-related effects on PCSK9 levels have a genetic basis. In addition, we compared the causal effects of PCSK9 on LDL-C levels between men and women and found a different response to statin treatment. This highlights the need for sex-sensitive dosages of lipid-lowering medication.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements.

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Characterization of Commercially Available Human Primary Alveolar Epithelial Cells.

In vitro lung research requires appropriate cell culture models that adequately mimic in vivo structure and function. Previously, researchers extensively used commercially available and easily expandable A549 and NCI-H441 cells, which replicate some but not all features of alveolar epithelial cells. Specifically, these cells are often restricted by terminally altered expression while lacking important alveolar epithelial characteristics. Of late, human primary alveolar epithelial cells (hPAEpCs) have become commercially available but are so far poorly specified. Here, we applied a comprehensive set of technologies to characterize their morphology, surface marker expression, transcriptomic profile, and functional properties. At optimized seeding numbers of 7,500 cells per square centimeter and growth at a gas-liquid interface, hPAEpCs formed regular monolayers with tight junctions and amiloride-sensitive transepithelial ion transport. Electron microscopy revealed lamellar body and microvilli formation characteristic for alveolar type II cells. Protein and single-cell transcriptomic analyses revealed expression of alveolar type I and type II cell markers; yet, transcriptomic data failed to detect NKX2-1, an important transcriptional regulator of alveolar cell differentiation. With increasing passage number, hPAEpCs transdifferentiated toward alveolar-basal intermediates characterized as SFTPC-, KRT8high, and KRT5- cells. In spite of marked changes in the transcriptome as a function of passaging, Uniform Manifold Approximation and Projection plots did not reveal major shifts in cell clusters, and epithelial permeability was unaffected. The present work delineates optimized culture conditions, cellular characteristics, and functional properties of commercially available hPAEpCs. hPAEpCs may provide a useful model system for studies on drug delivery, barrier function, and transepithelial ion transport in vitro.

Epigenome-Wide Meta-analysis Reveals Associations Between Dietary Glycemic Index and Glycemic Load and DNA Methylation in Children and Adolescents of Different Body Sizes.

OBJECTIVE: Dietary glycemic index (GI) and glycemic load (GL) are associated with cardiometabolic health in children and adolescents, with potential distinct effects in people with increased BMI. DNA methylation (DNAm) may mediate these effects. Thus, we conducted meta-analyses of epigenome-wide association studies (EWAS) between dietary GI and GL and blood DNAm of children and adolescents. RESEARCH DESIGN AND METHODS: We calculated dietary GI and GL and performed EWAS in children and adolescents (age range: 4.5-17 years) from six cohorts (N = 1,187). We performed stratified analyses of participants with normal weight (n = 801) or overweight or obesity (n = 386). We performed look-ups for the identified cytosine-phosphate-guanine (CpG) sites (false discovery rate [FDR] <0.05) with tissue-specific gene expression of 832 blood and 223 subcutaneous adipose tissue samples from children and adolescents. RESULTS: Dietary GL was positively associated with DNAm of cg20274553 (FDR <0.05), annotated to WDR27. Several CpGs were identified in the normal-weight (GI: 85; GL: 17) and overweight or obese (GI: 136; GL: 298; FDR <0.05) strata, and none overlapped between strata. In participants with overweight or obesity, identified CpGs were related to RNA expression of genes associated with impaired metabolism (e.g., FRAT1, CSF3). CONCLUSIONS: We identified 537 associations between dietary GI and GL and blood DNAm, mainly in children and adolescents with overweight or obesity. High-GI and/or -GL diets may influence epigenetic gene regulation and thereby promote metabolic derangements in young people with increased BMI.

An atlas of genome-wide gene expression and metabolite associations and possible mediation effects towards body mass index.

Investigating the cross talk of different omics layers is crucial to understand molecular pathomechanisms of metabolic diseases like obesity. Here, we present a large-scale association meta-analysis of genome-wide whole blood and peripheral blood mononuclear cell (PBMC) gene expressions profiled with Illumina HT12v4 microarrays and metabolite measurements from dried blood spots (DBS) characterized by targeted liquid chromatography tandem mass spectrometry (LC-MS/MS) in three large German cohort studies with up to 7706 samples. We found 37,295 associations comprising 72 amino acids (AA) and acylcarnitine (AC) metabolites (including ratios) and 8579 transcripts. We applied this catalogue of associations to investigate the impact of associating transcript-metabolite pairs on body mass index (BMI) as an example metabolic trait. This is achieved by conducting a comprehensive mediation analysis considering metabolites as mediators of gene expression effects and vice versa. We discovered large mediation networks comprising 27,023 potential mediation effects within 20,507 transcript-metabolite pairs. Resulting networks of highly connected (hub) transcripts and metabolites were leveraged to gain mechanistic insights into metabolic signaling pathways. In conclusion, here, we present the largest available multi-omics integration of genome-wide transcriptome data and metabolite data of amino acid and fatty acid metabolism and further leverage these findings to characterize potential mediation effects towards BMI proposing candidate mechanisms of obesity and related metabolic diseases. KEY MESSAGES: Thousands of associations of 72 amino acid and acylcarnitine metabolites and 8579 genes expand the knowledge of metabolome-transcriptome associations. A mediation analysis of effects on body mass index revealed large mediation networks of thousands of obesity-related gene-metabolite pairs. Highly connected, potentially mediating hub genes and metabolites enabled insight into obesity and related metabolic disease pathomechanisms.

Prognostic Factors for Postoperative Bleeding Complications and Prolonged Intensive Care after Percutaneous Hepatic Chemosaturation Procedures with Melphalan.

Percutaneous hepatic melphalan perfusion (chemosaturation) in patients with liver metastases is known to be associated with procedure-related hemodynamic depression and coagulation impairment, which may cause bleeding complications and/or a prolonged intensive care unit length of stay (ICU LOS). We retrospectively analyzed possible predictive factors for bleeding complications and an ICU LOS > 1 d in a cohort of 31 patients undergoing 90 chemosaturation procedures. Using a multivariable mixed-model approach, we identified the amount of perioperative fluid volume (OR 12.0, 95% CI 2.3-60.0, p = 0.003) and protamine (OR 0.065, 95% CI 0.007-0.55, p = 0.012) to be associated with bleeding complications. Furthermore, the amount of perioperative fluid volume was associated with an ICU LOS > 1 d (OR 5.2, 95% CI 1.4-19.0, p = 0.011). Heparin dosage, melphalan dosage, extracorporeal circulation time, and noradrenaline dosage had no significant effects on outcomes. Protamine use was not associated with anaphylactic or thromboembolic complications. Despite the limited sample size, these results suggest a restrictive perioperative fluid regime to be beneficial, and support the use of protamine for heparin reversal after chemosaturation procedures. Further prospective randomized trials are needed to confirm these findings.

SARS-CoV-2 Vaccine Breakthrough Infections of Omicron and Delta Variants in Healthcare Workers.

Understanding SARS-CoV-2 breakthrough infections in vaccinated healthcare workers is of key importance in mitigating the effects of the COVID-19 pandemic in healthcare facilities. An observational prospective cohort study was conducted in vaccinated employees with acute SARS-CoV-2 infection between October 2021 and February 2022. Serological and molecular testing was performed to determine SARS-CoV-2 viral load, lineage, antibody levels, and neutralizing antibody titers. A total of 571 (9.7%) employees experienced SARS-CoV-2 breakthrough infections during the enrolment period, of which 81 were included. The majority (n = 79, 97.5%) were symptomatic and most (n = 75, 92.6%) showed Ct values < 30 in RT-PCR assays. Twenty-four (30%) remained PCR-positive for > 15 days. Neutralizing antibody titers were strongest for the wildtype, intermediate for Delta, and lowest for Omicron variants. Omicron infections occurred at higher anti-RBD-IgG serum levels (p = 0.00001) and showed a trend for higher viral loads (p = 0.14, median Ct difference 4.3, 95% CI [-2.5-10.5]). For both variants, viral loads were significantly higher in participants with lower anti-RBD-IgG serum levels (p = 0.02). In conclusion, while the clinical course of infection with both the Omicron and Delta variants was predominantly mild to moderate in our study population, waning immune response over time and prolonged viral shedding were observed.

Imputation-powered whole-exome analysis identifies genes associated with kidney function and disease in the UK Biobank.

Genome-wide association studies have discovered hundreds of associations between common genotypes and kidney function but cannot comprehensively investigate rare coding variants. Here, we apply a genotype imputation approach to whole exome sequencing data from the UK Biobank to increase sample size from 166,891 to 408,511. We detect 158 rare variants and 105 genes significantly associated with one or more of five kidney function traits, including genes not previously linked to kidney disease in humans. The imputation-powered findings derive support from clinical record-based kidney disease information, such as for a previously unreported splice allele in PKD2, and from functional studies of a previously unreported frameshift allele in CLDN10. This cost-efficient approach boosts statistical power to detect and characterize both known and novel disease susceptibility variants and genes, can be generalized to larger future studies, and generates a comprehensive resource ( https://ckdgen-ukbb.gm.eurac.edu/ ) to direct experimental and clinical studies of kidney disease.

Epidemiological Analysis of the Emergency Vascular Access in Pediatric Trauma Patients: Single-Center Experience of Intravenous, Intraosseous, Central Venous, and Arterial Line Placements.

Vascular access in severely injured pediatric trauma patients is associated with time-critical circumstances and low incidences, whereas only scarce literature on procedure performance is available. The purpose of this study was to analyze the performance of different vascular access procedures from the first contact at the scene until three hours after admission. Intubated pediatric trauma patients admitted from the scene to a single Level I trauma center between 2008 and 2019 were analyzed regarding intravenous (IV) and intraosseous (IO) accesses, central venous catheterization (CVC) and arterial line placement. Sixty-five children with a median age of 14 years and median injury severity score of 29 points were included, of which 62 (96.6%) underwent successful prehospital IV or IO access by emergency medical service (EMS) physicians, while it failed in two children (3.1%). On emergency department (ED) admission, IV cannulas of prehospital EMS had malfunctions or were dislodged in seven of 55 children (12.7%). IO access was performed in 17 children without complications, and was associated with younger age, higher injury severity and higher mortality. Fifty-two CVC placements (58 attempts) and 55 arterial line placements (59 attempts) were performed in 45 and 52 children, respectively. All CVC and arterial line placements were performed in the ED, operating room (OR) and intensive care unit (ICU). Ten mechanical complications related to CVC placement (17.8%) and seven related to arterial line placement (10.2%) were observed, none of which had outcome-relevant consequences. This case series suggests that mechanical issues of vascular access may frequently occur, underlining the need for special preparedness in prehospital, ED, ICU and OR environments.

MicroRNA-223 Dampens Pulmonary Inflammation during Pneumococcal Pneumonia.

Community-acquired pneumonia remains a major contributor to global communicable disease-mediated mortality. Neutrophils play a leading role in trying to contain bacterial lung infection, but they also drive detrimental pulmonary inflammation, when dysregulated. Here we aimed at understanding the role of microRNA-223 in orchestrating pulmonary inflammation during pneumococcal pneumonia. Serum microRNA-223 was measured in patients with pneumococcal pneumonia and in healthy subjects. Pulmonary inflammation in wild-type and microRNA-223-knockout mice was assessed in terms of disease course, histopathology, cellular recruitment and evaluation of inflammatory protein and gene signatures following pneumococcal infection. Low levels of serum microRNA-223 correlated with increased disease severity in pneumococcal pneumonia patients. Prolonged neutrophilic influx into the lungs and alveolar spaces was detected in pneumococci-infected microRNA-223-knockout mice, possibly accounting for aggravated histopathology and acute lung injury. Expression of microRNA-223 in wild-type mice was induced by pneumococcal infection in a time-dependent manner in whole lungs and lung neutrophils. Single-cell transcriptome analyses of murine lungs revealed a unique profile of antimicrobial and cellular maturation genes that are dysregulated in neutrophils lacking microRNA-223. Taken together, low levels of microRNA-223 in human pneumonia patient serum were associated with increased disease severity, whilst its absence provoked dysregulation of the neutrophil transcriptome in murine pneumococcal pneumonia.

National and subnational short-term forecasting of COVID-19 in Germany and Poland during early 2021.

BACKGROUND: During the COVID-19 pandemic there has been a strong interest in forecasts of the short-term development of epidemiological indicators to inform decision makers. In this study we evaluate probabilistic real-time predictions of confirmed cases and deaths from COVID-19 in Germany and Poland for the period from January through April 2021. METHODS: We evaluate probabilistic real-time predictions of confirmed cases and deaths from COVID-19 in Germany and Poland. These were issued by 15 different forecasting models, run by independent research teams. Moreover, we study the performance of combined ensemble forecasts. Evaluation of probabilistic forecasts is based on proper scoring rules, along with interval coverage proportions to assess calibration. The presented work is part of a pre-registered evaluation study. RESULTS: We find that many, though not all, models outperform a simple baseline model up to four weeks ahead for the considered targets. Ensemble methods show very good relative performance. The addressed time period is characterized by rather stable non-pharmaceutical interventions in both countries, making short-term predictions more straightforward than in previous periods. However, major trend changes in reported cases, like the rebound in cases due to the rise of the B.1.1.7 (Alpha) variant in March 2021, prove challenging to predict. CONCLUSIONS: Multi-model approaches can help to improve the performance of epidemiological forecasts. However, while death numbers can be predicted with some success based on current case and hospitalization data, predictability of case numbers remains low beyond quite short time horizons. Additional data sources including sequencing and mobility data, which were not extensively used in the present study, may help to improve performance.

We compare forecasts of weekly case and death numbers for COVID-19 in Germany and Poland based on 15 different modelling approaches. These cover the period from January to April 2021 and address numbers of cases and deaths one and two weeks into the future, along with the respective uncertainties. We find that combining different forecasts into one forecast can enable better predictions. However, case numbers over longer periods were challenging to predict. Additional data sources, such as information about different versions of the SARS-CoV-2 virus present in the population, might improve forecasts in the future.

On the Parametrization of Epidemiologic Models-Lessons from Modelling COVID-19 Epidemic.

Numerous prediction models of SARS-CoV-2 pandemic were proposed in the past. Unknown parameters of these models are often estimated based on observational data. However, lag in case-reporting, changing testing policy or incompleteness of data lead to biased estimates. Moreover, parametrization is time-dependent due to changing age-structures, emerging virus variants, non-pharmaceutical interventions, and vaccination programs. To cover these aspects, we propose a principled approach to parametrize a SIR-type epidemiologic model by embedding it as a hidden layer into an input-output non-linear dynamical system (IO-NLDS). Observable data are coupled to hidden states of the model by appropriate data models considering possible biases of the data. This includes data issues such as known delays or biases in reporting. We estimate model parameters including their time-dependence by a Bayesian knowledge synthesis process considering parameter ranges derived from external studies as prior information. We applied this approach on a specific SIR-type model and data of Germany and Saxony demonstrating good prediction performances. Our approach can estimate and compare the relative effectiveness of non-pharmaceutical interventions and provide scenarios of the future course of the epidemic under specified conditions. It can be translated to other data sets, i.e., other countries and other SIR-type models.

Chest Tube Placement in Mechanically Ventilated Trauma Patients: Differences between Computed Tomography-Based Indication and Clinical Decision.

The rate of occult pneumothorax in intubated and mechanically ventilated trauma patients until initial computed tomography (CT) remains undetermined. The primary aims of this study were to analyze initial chest CTs with respect to the thoracic pathology of trauma, the clinical injury severity, and chest tube placement (CTP) before and after CT. In a single-center retrospective analysis of 616 intubated and mechanically ventilated adult patients admitted directly from the scene to the emergency department (ED), 224 underwent CTP (36%). Of these, 142 patients (62%) underwent CTP before CT, of which, 125 (88%) had significant chest injury on CT. Seventeen patients had minor or absent chest injuries, most of which were associated with transient or unrecognized tracheal tube malposition. After CT, CTP was performed in another 82 patients, of which, 56 (68.3%) had relevant pneumothorax and 26 had minor findings on CT. Sixty patients who had already undergone CTP before CT received another CTP after CT, of which, 15 (25%) had relevant pneumothorax and 45 (75%) had functionality issues or malposition requiring replacement. Nine patients showed small pneumothorax on CT, and did not undergo CTP (including four patients with CTP before CT). The physiological variables were unspecific, and the trauma scores were dependent on the CT findings for identifying patients at risk for CTP. In conclusion, the clinical decisions for CTP before CT are associated with relevant false-negative and false-positive cases. Clinical assessment and CT imaging, together, are important indicators for CTP decisions that cannot be achieved by using clinical assessment or CT alone.

First genome-wide association study of 99 body measures derived from 3-dimensional body scans.

Body height, body mass index, hip and waist circumference are important risk factors or outcome variables in clinical and epidemiological research with complex underlying genetics. However, these classical anthropometric traits represent only a very limited view on the human body and other traits with potentially higher functional specificity are not yet studied to a larger extent. Participants of LIFE-Adult were assessed by three-dimensional body scanner VITUS XXL determining 99 high-quality anthropometric traits in parallel. Genotyping was performed by Axiom Genome-Wide CEU 1 Array Plate microarray technology and imputation was done using 1000 Genomes phase 3 reference panel. Combined phenotype and genetic information are available for a total of 7,562 participants. Largest heritabilities were estimated for height traits (maximum heritability with h2 = 44% for neck height) and 61 traits achieved values larger than 20%. By genome-wide analyses, we identified 16 loci associated with at least one of the 99 traits. Ten of these loci were not described for association with classical anthropometric traits so far. The strongest novel association was observed for 7p14.3 (rs11979006, P = 2.12 × 10-9) for the trait Back Width with ZNRF2 as the most plausible candidate gene. Loci established for association with classical anthropometric traits were subjected to anthropometric phenome-wide association analysis. From the reported 709 loci, 211 are co-associated with body scanner traits (enrichment: OR = 1.96, P = 1.08 × 10-61). We conclude that genetics of 3D laser-based anthropometry is promising to identify novel loci and to improve the functional understanding of established ones.

A pulmonologist's guide to perform and analyse cross-species single lung cell transcriptomics.

Single-cell ribonucleic acid sequencing is becoming widely employed to study biological processes at a novel resolution depth. The ability to analyse transcriptomes of multiple heterogeneous cell types in parallel is especially valuable for cell-focused lung research where a variety of resident and recruited cells are essential for maintaining organ functionality. We compared the single-cell transcriptomes from publicly available and unpublished datasets of the lungs in six different species: human (Homo sapiens), African green monkey (Chlorocebus sabaeus), pig (Sus domesticus), hamster (Mesocricetus auratus), rat (Rattus norvegicus) and mouse (Mus musculus) by employing RNA velocity and intercellular communication based on ligand-receptor co-expression, among other techniques. Specifically, we demonstrated a workflow for interspecies data integration, applied a single unified gene nomenclature, performed cell-specific clustering and identified marker genes for each species. Overall, integrative approaches combining newly sequenced as well as publicly available datasets could help identify species-specific transcriptomic signatures in both healthy and diseased lung tissue and select appropriate models for future respiratory research.

Relationship between impaired BMP signalling and clinical risk factors at early-stage vascular injury in the preterm infant.

INTRODUCTION: Chronic lung disease, that is, bronchopulmonary dysplasia (BPD) is the most common complication in preterm infants and develops as a consequence of the misguided formation of the gas-exchange area undergoing prenatal and postnatal injury. Subsequent vascular disease and its progression into pulmonary arterial hypertension critically determines long-term outcome in the BPD infant but lacks identification of early, disease-defining changes. METHODS: We link impaired bone morphogenetic protein (BMP) signalling to the earliest onset of vascular pathology in the human preterm lung and delineate the specific effects of the most prevalent prenatal and postnatal clinical risk factors for lung injury mimicking clinically relevant conditions in a multilayered animal model using wild-type and transgenic neonatal mice. RESULTS: We demonstrate (1) the significant reduction in BMP receptor 2 (BMPR2) expression at the onset of vascular pathology in the lung of preterm infants, later mirrored by reduced plasma BMP protein levels in infants with developing BPD, (2) the rapid impairment (and persistent change) of BMPR2 signalling on postnatal exposure to hyperoxia and mechanical ventilation, aggravated by prenatal cigarette smoke in a preclinical mouse model and (3) a link to defective alveolar septation and matrix remodelling through platelet derived growth factor-receptor alpha deficiency. In a treatment approach, we partially reversed vascular pathology by BMPR2-targeted treatment with FK506 in vitro and in vivo. CONCLUSION: We identified impaired BMP signalling as a hallmark of early vascular disease in the injured neonatal lung while outlining its promising potential as a future biomarker or therapeutic target in this growing, high-risk patient population.

Genetically regulated gene expression and proteins revealed discordant effects.

BACKGROUND: Although gene-expression (GE) and protein levels are typically strongly genetically regulated, their correlation is known to be low. Here we investigate this phenomenon by focusing on the genetic background of this correlation in order to understand the similarities and differences in the genetic regulation of these omics layers. METHODS AND RESULTS: We performed locus-wide association studies of 92 protein levels measured in whole blood for 2,014 samples of European ancestry and found that 66 are genetically regulated. Three female- and one male-specific effects were detected. We estimated the genetically regulated GE for all significant genes in 49 GTEx v8 tissues. A total of 7 proteins showed negative correlations with their respective GE across multiple tissues. Finally, we tested for causal links of GE on protein expression via Mendelian Randomization, and confirmed a negative causal effect of GE on protein level for five of these genes in a total of 63 gene-tissue pairs: BLMH, CASP3, CXCL16, IL6R, and SFTPD. For IL6R, we replicated the negative causal effect on coronary-artery disease (CAD), while its GE was positively linked to CAD. CONCLUSION: While total GE and protein levels are only weakly correlated, we found high correlations between their genetically regulated components across multiple tissues. Of note, strong negative causal effects of tissue-specific GE on five protein levels were detected. Causal network analyses revealed that GE effects on CAD risks was in general mediated by protein levels.

Transcriptome Analyses of Adipose Tissue Samples Identify EGFL6 as a Candidate Gene Involved in Obesity-Related Adipose Tissue Dysfunction in Children.

Obesity develops early in childhood and is accompanied by early signs of adipose tissue (AT) dysfunction and metabolic disease in children. In order to analyse the molecular processes during obesity-related AT accumulation in children, we investigated genome-wide expression profiles in AT samples, isolated adipocytes, and stromal vascular fraction (SVF) cells and assessed their relation to obesity as well as biological and functional AT parameters. We detected alterations in gene expression associated with obesity and related parameters, i.e., BMI SDS, adipocyte size, macrophage infiltration, adiponectin, and/or leptin. While differential gene expression in AT and adipocytes shared an enrichment in metabolic pathways and pathways related to extracellular structural organisation, SVF cells showed an overrepresentation in inflammatory pathways. In adipocytes, we found the strongest positive association for epidermal growth factor-like protein 6 (EGFL6) with adipocyte hypertrophy. EGFL6 was also upregulated during in vitro adipocyte differentiation. In children, EGFL6 expression was positively correlated to parameters of AT dysfunction and metabolic disease such as macrophage infiltration into AT, hs-CRP, leptin levels, and HOMA-IR. In conclusion, we provide evidence for early alterations in AT gene expression related to AT dysfunction in children and identified EGFL6 as potentially being involved in processes underlying the pathogenesis of metabolic disease.

Impact of medication on blood transcriptome reveals off-target regulations of beta-blockers.

BACKGROUND: For many drugs, mechanisms of action with regard to desired effects and/or unwanted side effects are only incompletely understood. To investigate possible pleiotropic effects and respective molecular mechanisms, we describe here a catalogue of commonly used drugs and their impact on the blood transcriptome. METHODS AND RESULTS: From a population-based cohort in Germany (LIFE-Adult), we collected genome-wide gene-expression data in whole blood using in Illumina HT12v4 micro-arrays (n = 3,378; 19,974 gene expression probes per individual). Expression profiles were correlated with the intake of active substances as assessed by participants' medication. This resulted in a catalogue of fourteen substances that were identified as associated with differential gene expression for a total of 534 genes. As an independent replication cohort, an observational study of patients with suspected or confirmed stable coronary artery disease (CAD) or myocardial infarction (LIFE-Heart, n = 3,008, 19,966 gene expression probes per individual) was employed. Notably, we were able to replicate differential gene expression for three active substances affecting 80 genes in peripheral blood mononuclear cells (carvedilol: 25; prednisolone: 17; timolol: 38). Additionally, using gene ontology enrichment analysis, we demonstrated for timolol a significant enrichment in 23 pathways, 19 of them including either GPER1 or PDE4B. In the case of carvedilol, we showed that, beside genes with well-established association with hypertension (GPER1, PDE4B and TNFAIP3), the drug also affects genes that are only indirectly linked to hypertension due to their effects on artery walls or their role in lipid biosynthesis. CONCLUSIONS: Our developed catalogue of blood gene expressions profiles affected by medication can be used to support both, drug repurposing and the identification of possible off-target effects.